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Abstract 

The nucleon-nucleon scattering in a large magnetic background is considered 
to find its potential to change the neutrino emissivity of the neutron stars. For 
this purpose we consider the one-pion-exchange approximation to find the NN 
cross-section in a background field as large as 10 15 G — 10 18 G. We show that 
the NN cross-section in neutron stars with temperatures in the range 0.1-5 MeV 
can be changed up to the one order of magnitude with respect to the one in 
the absence of the magnetic field. In the limit of the soft neutrino emission 
the neutrino emissivity can be written in terms of the NN scattering amplitude 
therefore the large magnetic fields can dramatically change the neutrino emissivity 
of the neutron stars as well. 
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1 Introduction 



Neutron stars in general can have magnetic fields as large as 10 12 G — 10 18 G [T], [2], [3], 
II], 0, [6] while the timing observations indicate that the surface magnetic fields of the 
radio pulsars are about 10 12 G [7], by equating magnetic field energy to the gravitational 
binding energy, one can estimate an upper limit for strength of the magnetic field to 
be B ~ 10 18 G [8]. Meanwhile cooling processes in the neutron stars formation that 
has long been considered by many physicists, can be affected by such a large magnetic 
field in two ways. It can change the electron or the photon conductivities [9], [2], [10] as 
well as the wave functions and the propagators of the particles involving in the cooling 
processes. One of the most important reactions during the thermal evolution of the 
neutron stars is the emission of neutrino via the bremsstrahlung process 

NN — >■ NN + vu. (1) 

In the first sight it seems that since the mass of the nucleons are about lGev therefore 
the magnetic fields even as large as 10 18 G can not significantly change the cross sec- 
tion. However one should note that there are more energy scales besides the mass of 
the nucleons. In fact one of the scales is the temperature of the evolving neutron star 
which becomes less than lMeV in a few seconds of the late-time cooling phase and is 
comparable with B ~ 10 13 G. Furthermore if one consider one-pion-exchange (OPE) for 
the nucleon-nucleon interaction then the magnetic fields of the order of 10 17 G and more 
can change the strong coupling constant as well. Meanwhile in such a low temperature 
era the bremsstrahlung process can only produce the soft neutrinos with respect to 
the initial energies and therefore the neutrino emissivity can be related to the on-shell 
NN scattering amplitude [IT] . Thus to find the neutrino emissivity one can study the 
nucleon-nucleon scattering separately. The amplitude can be even constructed from 
NN scattering data, as is done in [11] as a model independent way to describe the 
neutrino emissivity. But the usual studies on the nucleon-nucleon scattering in OPE 
approximation |12j and more precisely the existing data are in the absence of the mag- 
netic field. Our aim in this paper is to examine the magnetic field dependence of the 
nucleon-nucleon scattering in the OPE approximation when soft neutrino approxima- 
tion is applicable. 

The paper is organized as follows: in section 2 we give the field operators of the 
proton and neutrons in an external magnetic field. Consequently we obtain the cross 
section of the scattering in the magnetic field in section 3. An estimation for the effect 
of strong magnetic field on the neutrino emissivity is presented in section 4. In section 5 
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we give some concluding remarks. Appendix A is devoted to solving the Dirac equation 
for proton and neutron in an external magnetic field. In Appendix B we present a 
simplified expression for the amplitude matrix elements of the scattering process. 



2 Fermion field operator 

To study the proton-neutron scattering in a strong external magnetic field, the field 
operators of the proton and neutron are needed. Since the field operators contain wave 
functions and energy spectrum, Dirac equation in presence of an external magnetic 
field must be solved. We choose the external magnetic field to be in the positive z-axis 
B = Bk. Because of gauge invariance there is a freedom in the choice of vector potential 
Apix), thus for simplicity we choose 

A^x) = (0,yB,0,0)- (2) 

This choice of vector potential pushes all the ^/-dependence into the spinors alone and 
not in the phases, so it is useful to introduce the following notations 

x i = (t,x,Q,z), 

% = (V X ,V S ), (3) 

where is any 2-vector. The solutions of the Dirac equation for proton and neutron in 
the external magnetic field (j2J) are presented in the Appendix A. By using this solutions 
the proton field operator can be written as 

n=0 s=±l ^ ^ ' 

(4) 

where the creation and annihilation operators obey the following anticommutation re- 
lations 

{a s (n,pj),al,(n',^)} = {b s (n,p^),bl,(n',^)} = {2^) 2 5 s ^5 n)nl 5 2 i {p - p). (5) 

With the help of (JSJ), the explicit form of the wave functions given in the equations 
([ST]) . (IMjl and the completeness relation (155)1 . it can be shown that the equal-time 
anticommutation relation is satisfied 

U; p (t,x)^l(t,x')} = 5 3 (x-x'). (6) 
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For convenience we normalize our 1-particle states in a box with dimensions L x L y L z = V 
such that the proton state is 

\p{n,pj,s)>= J— al(n,p^ | > . (7) 

V ^x^z 

Similarly for the neutron field operator one has 

M*) = E / ^N s [aMUs(p)e^ x + bl{p)V s {f)e +i n, (8) 

in which the creation and annihilation operators obey the following anticommutation 
relations 

= {b s (p),bl@)} = (2tt)V*V -£)■ (9) 
By using the relations and the explicit form of the wave functions given in the 
equations fl66l) and fl69|) it is easy to show that equal-time anticommutation relation is 
satisfied 

{Mt,x)^i(t,x')} = 5 3 (x-x'). (10) 
Box normalization condition for the neutron 1-particle state is 

\n{p,s)>=^=at(p) |0>. (11) 



3 Cross section of the scattering in the external 
magnetic field 

Our aim in this section is to explore the scattering cross section of p + n — > p + n in the 
background of a strong magnetic field, erg, using the one-pion-exchange approximation. 
The deferential cross section can be written as 



V I S 1 ! 2 

da = p ^r l -+dp } (12) 

\v„ - v n T 



v 



where \v p — v n \/V is the flux of the incident particles, T is the time interval, \S\ 2 is the 
average of the squared scattering matrix elements and dp is the deferential phase space 
of the final states. The effective low energy Lagrangian for pion-nucleon interaction is 

£int = g^N^pl^Ipp^ + V / 2'0 n 7 5 '0p0 7 r- + y/2^ p ^1p n (f)^ - ^nj^n&n ], (13) 

where — 14 is the pion-nucleon coupling constant and ip n , ip p and are neutron, 
proton and pion fields, respectively. To obtain the scattering matrix elements in the 
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external magnetic field we use the standard method where the wave functions and 
propagators are modified with the external magnetic field while the vertices are leaved 
unchanged. The proton-neutron scattering at the tree level involves two diagrams, one 
is mediated by the neutral pion S nc 



S nc =< p'{n',^,s'),n'(k',r')\T{ J d x'd xip p (x')(ig nN j )ifj p (x')<j) n o(x') (14) 

4>n(x)(-ignNl 5 )*Pn(x)(pKo(x)}\p(n,p^ s),7l(k,r) >, 

and the other is mediated by the charged pion S cc 

S cc =<p'(n'J^s'),n'(k',r')\T{ j rfVrf 4 x^(x , )(^v / 2^ J v7 5 )^n(x , )0. + (x , ) (15) 

^n(x)(iV2g wN j 5 )i! p (x)(p 7T -(x)}\p(n,p^,s),n(k,r) >, 

where p (n) and p' (n'), respectively, correspond to the incoming and outgoing protons 
(neutrons). Since the neutral pion has no charge its propagator in the external magnetic 
field does not change, but the propagator of charged pion will be modified. In fact the 
propagator of the charged pion as a charged scalar particle in an external magnetic field 
can be obtained as 1 13 1 



f*x "i f ,j4 

< T{^ + (x')^-(x)} >= expHe J d£p{A" + -F""(£ - *)„}) J -X e -*«V-«) 

X / t S R \ ex P(^{^ll - tan( R jBg) gl - m '})>(16) 
J cos{eBs) " eBs 

in which is the field strength tensor of the external magnetic field and 

4 = (qo) 2 - (g 3 ) 2 , ql = {qif + fe) 2 . (17) 

The temperature of a neutron star in a few seconds in the late-time cooling phase, 
drops to less than 1 MeV therefore the energy is sufficiently low to use the four-fermion 
approximation, q 2 /m% = 0. In this case the propagator (|T6|) can be cast into 

<T{ ( f> 7T+ (x')(t> 7r -(x)}> = ^D(eB/ml)5 A (x' - x); 

7T 

poo -s 

D(eB/ml) = / ds (18) 
J cosh(^) 

In the four-fermion approximation, D actually shows how the strong coupling constant 
depends on the external magnetic field. The B-dependence of the coupling constant 
is explicitly shown in Fig. 1. By using the explicit form of the proton and neutron 
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Figure 1: The B- dependence of the squared coupling constant. The variation starts 
from 10 15 G but the value decreases significantly from 10 17 G. 



field operators (TjJ and ([8]) and the definitions of the 1-particle states ([7]) and ffTTj) . the 
scattering matrix elements can be obtained as follows 



nc(cc) 



^ p ^ n ^ p /^ n /(27^) J ^J(p , + k' — p — k)(-iM nc{cc) ), 



L X L Z V 

where the neutral current scattering amplitude is given by 



(19) 



M, 



n 2 



nc n 

mi 



dye^K-^v [U p , (y, n', pj, s')j%(y, n, p+ a)) [U n , (k', r'^U^k, r)} , (20) 



and the charged current is 



Mr, = 2D 



n 2 
m 2 



dye-^ k 'y- k y^[U p ,(y, n'J p s'h 5 U n (k, r)][U n ,(k\ r')j%(y, n,p^, s)].(21) 



Now we are ready to calculate the scattering cross section in the external magnetic field. 
For this purpose we have 
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(2s 



> + l)(2s n + l)^ lSnc + S ' 

' v ' spin 



2 

cc | ) 



(22) 



where using the normalization relation for the Dirac 5-functions 

1 



3/' 



p' + k' — p — k)Y 



(2tt) 3 
6 



TL x L z 8l(p' + k'-p-k), 



(23) 



cast the squared amplitude into 

\S\ 2 = ^^WW(27r) 8 < §(p' + A/-p-A ; )|A<| 2 . (24) 

Now the equations (fT2|) . (124|) and the deferential phase space of the final states in the 
external magnetic field 

d 2 v'i d 3 k' 

^Wf-^W (25) 

can be used to find the corresponding deferential cross section as follows 

B |^-4|(2 Sp + l)(2 Sn + l) s ^ p n (27r) 2 (2vr)3 »' 

x(2n) 3 5l(p' + k' -p- k)\M\ 2 . (26) 

Since the nucleons in the neutron star are nonrelativistic their energies and wave func- 
tions can be expanded up to the order of 0(l/mjv) and the amplitude M. can be 
obtained analytically as is given in the Appendix B. Meanwhile one can easily show 
that 

N p , N n , N p r, N n , = l + 0(l/m 2 N ). (27) 
Therefore the leading term in the scattering cross section (|26|) is 

d% d 3 k' 



_ it I 



, 2! , +1 )' (2! , +1) £My^->-»i^i ! (* 



(2s p + l)(2s n + l) 

v ^ * v ' spins 

After integrating over p^., fc^ with the help of the Dirac ^-functions and defining 

k' x = k' ± cos &, k' y = k' ± sin &, (29) 

equation fT28|) leads to 

I 00 r+00 p2tt r+00 

° b = ^ K -v n \^ o L dp '-L M 'l dk ' A 
x (30) 

v ^ ' v ' spins 

By using the nonrelativistic expansion of the energies as are given in the equations (1761) 
and (1771) one has 

= E' p + E' n — E p — E n 

k'l + 2(p' z -l( Pz + k z )f - 2eB(n' max - n')) , (31) 



2m N V 2 



where n' max is defined as 

n' max ^n+^(s'-s) + ^(r'-r) + ^(kl + \{p z - k z f). (32) 
The positive solution of (I3"TT) is 



±0 = y 2eB(n' max - nf) - 2(p' z - l{p M + k z )) 2 . (33) 
Now one can perform the integration over k' ± to find 

" = 4^J(5^WT) E E /. **/ < 34 > 

^ ' v ' spins ri—0 z- 

where 

A^^INTKJ, (35) 
while the integration bounds over p' z can be determined by the following equation 

K± = l( Pz + k z )± ^eB{n' max -n'). (36) 



To evaluate Ob in the equation (|34l) first one needs to determine the maximum Landau- 
level that proton can occupy in the initial state, i.e. M ma x- Since we would like to 
consider the process (JT|) in the temperature range T ~ lMev [2], [I], [TJ] therefore 
nucleons are nonrelativistic and the energy eigenvalues for the protons and neutrons 
can be given by ( l76p and fl77j) . respectively. According to the Virial theorem, each 
nucleon at a temperature T has a momentum for each component about (mjyT) 1 / 2 . 
Therefore one can easily obtain the maximum Landau-level by comparing the magnetic 
energy contribution to the proton neB/rriN, see (J76l) . with the thermal energy as N max = 
INT(m N T/eB). 

In this section we obtained the cross-section in the background of a strong magnetic 
field erg by considering OPE in the four-fermion approximation. For simplicity we 
calculate the cross-section in the head-on collisions i.e. (p + k)j = 0. Since we would 
like to find if the magnetic field can significantly change the thermal evolution of the 
neutron stars, it would be better to normalize our result with the cross-section in the 
absence of the magnetic field, <To, but with the same approximations and conditions i.e. 
(p + k)= 0, or 

?>Timl\v p -v n \m% 
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2.5 




Figure 2: The variation of the T in terms of Log(B/lG) for different temperatures. 
Therefore the deviation — as a function of the magnetic field and temperature can be 



given by 



-i max q -i A/" max 

V(B T\ = V ° B = ( mN \5 1 V 

1 ' ' ' (l+^ma X ) ^ CTO 87T 1 \p\ } (1 + Nmax) ^ 

x x x 



'"max Pr , J / /*Z7T 



(2s p + l)(2s n + 1) 4^ ^ Jp' 7 

^ v spins n—0 ' 



\M\* = -±\M\*. (39) 



where 

~|2 - A/1 ! - 

The variation of T as a function of the external magnetic field for different temperatures 
are shown in the Figures 2 and 3. 



4 Neutrino emissivity 

The NN-scattering cross section in the magnetic fields as large as 10 15 G — 10 18 G shows 
that the strong magnetic filed of the neutron stars might change the neutrino emissivity 
and therefore the cooling process. The rate of emission and absorbtion of neutrino in the 
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Figure 3: The variation of the F in terms of Log(B/lG) for different temperatures. 

strong magnetic field of a neutron star have been already considered for the Urea process 
in [15], [16J and [IT]. In [15] for a strange stars with all charged particles in the lowest 
Landau level, it is shown that the magnetic field accelerates the cooling process. While 
in [16] it is found that the direct Urea emissivity is strongly enhanced and oscillates as 
a function of the magnetic field which globally tends to average out. Here we would like 
to consider the effect of strong magnetic field on the neutrino bremsstrahlung process 
([1]). To obtain the emissivity at the tree level one should consider all Feynman diagrams 
at this level. Here for simplicity we just consider the diagram given in FigJH to give 
an order of magnitude for the emissivity and leave the complete calculation for future 
study. For this purpose we use the pion-nucleon interaction Lagrangian fflBl and the 
weak interactions Lagrangian 

Ant = -^[^7 M (l-7 5 )^-^(l-7 5 )^n + ^(l-4sin 2 ^- 7 5 )^]^. (40) 
4 cos & w p 

According to the Fermi's golden rule, neutrino emissivity in the magnetic field can be 
written as 

^ oo^ oo^ d 2 p ^ LxLz d 3 ky d 2 p'^L x L z d^k'V d 3 qi V d 3 q 2 V 
eB ~VThhJ (2vr) 2 (2tt)3 (2vrP (2vf (2tt)3 (2^ 

where u is the total energy of neutrino with the momentum of qi and antineutrino with 
the momentum of q 2 , all other momenta have the same definition as was given in the 



5>i 



u, (41) 



\sptns 
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Figure 4: Feynman diagram for neutrino emissivity. 



last section. F is defined in terms of the Fermi-Dirac distribution functions of nucleons 

as 

^ = /p/»(l-/pO(l-/nO; /=: 1 - 7 e^, (42) 

1 + exp(^^) 

in which \x is the chemical potential of the nucleon. Now we consider the nucleons, 
comparing the available energies, as non-relativistic particles. Furthermore the anoma- 
lous magnetic moment in the energy spectrum of nucleons are negligible, see Eqs. (1761) . 
(j77|) . Meanwhile since the neutrinos are soft the three- momenta of neutrino and an- 
tineutrino can be dropped from the Dirac 5-functions. Therefore the spin summed over 
the squared scattering matrix element can be obtained as 

£ |S| ^ (U^k (2 " ! *' + k' + q,+ q 2-p-k)Y 1 \M?, (43) 

spin v 1 spin 

where 

E \ M I 2 = (^rK) 2 7^TFyM(p' z - P>? + 2eB(n' + n)) 
f-* n \V2m 2 N mlJ (E 1 + E 2 ) 2 

x(-fn,n'-C„' + In-l,n>-ll* n -l,n<-l) ~ ^BVr}fn{In,n'I* n -l,n'-l + 4-l,n'-l^„/)}, (44) 

in which the functions l n . m are defined in fTTTT) . Substituting (|43p in the emissivity 
expression (T4"T1) yields 

X X ' d 2 Vi d 3 k d 2 ^ <Pk> d 3 Ql 1 d 3 q 2 1 



6 B = v^bJ2J2 

n=0 n'=0 



(2vr) 2 (2tt) 3 (2tt) 2 (2vr) 3 (2tt) 3 2^ (2vr) 3 2E 2 



x(2n) 3 5l(p' + k' + gi + q 2 - p - k)F I \M\Au:. ( 15) 
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T=3Mev 
T=5Mev 
T=7Mev 



17.6 17.8 18 18.2 18.4 18.6 18.8 19 19.2 
Log(B/lG) 



Figure 5: Y = ^ is plotted as function of magnetic field for different temperatures. 

One should note that the range of L y can be fixed in a natural way as L y = ^= 
[15] . Furthermore to perform the integrals the momentum of particles for a degenerate 
system, can not be larger than the Fermi momenta that in turns depend on the density 
of system and implicitly on the temperature [12] . so the phase space integrals of the 
initial particles confined to the Fermi surface. While the remanning integrals on the 
phase space of the final particles can be accomplished using the momentum conserving 
Dirac 5-functions . Here we define T = ^ where e shows the corresponding emissivity 
in the absence of magnetic field. The variation of T with respect to the neutron star 
magnetic field for different temperature is shown in Figj5l 

5 Conclusion 

We have examined the nucleon-nucleon scattering in a large magnetic background in a 
neutron star at the temperature about 1 MeV . Since the energy is low we can safely use 
the one-pion-exchange approximation and the four-fermion approximation. In fact the 
external magnetic field can change the effective coupling constant that is shown in the 
Fig.l. As the figure shows the coupling is insensitive to the background field up to the 
10 17 G and then by increasing the magnetic field to 4 x 10 18 G the squared value of the 
strong coupling constant decreases to 0.6 its initial value. However the magnetic field 
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can change the nucleon wave functions and their corresponding field operators as well, 
see equations (4) and (8). Therefore the cross-section has an overall change as given in 
equation (34). The variation of the cross-section, for head-on collisions, in a magnetic 
background that is normalized with its correspondence without background but with 
the same approximations are drawn in the Fig. 2 and Fig. 3 for different temperatures. 
As the figures show the deviation for the lower temperature start in the lower magnetic 
fields. Nevertheless to have a significant deviation the magnetic field even for the 
temperature as low as 0.1 MeV, should be grater than 10 15 G. We also studied the neutrino 
bremsstrahlung emissivity by considering just one Feynman diagram for this process at 
the tree level to give an estimate on the neutrino emissivity in a strong magnetic field. 
As is shown in the Figj5] around B = 10 18 G the emissivity changes about an order of 
magnitude with respect to its corresponding value without magnetic filed. Furthermore 
the enhancement on the emissivity is higher for the lower temperature. Therefore the 
thermal evolution of the neutron stars in the late time cooling era for B > 10 15 G can 
be dramatically changed. 
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A Solutions of Dirac equation 

In this Appendix we drive the wave functions of proton and neutron in the external 
magnetic field by tracing a procedure similar to [18] . Dirac equation for a particle with 
mass m, charge e = \e\ and anomalous magnetic moment pL = K-£- in presence of the 
external magnetic field (EJ) is 



dip 

Where the canonical momentum II is given by the minimal coupling 

II = -zV - eA. 

By using the standard Pauli-Dirac representation for Dirac matrices 



(46) 



(47) 



a 



a 
a 



I 
—I 



<T Z 

o z 



I 

1 



, (4* 



-ip.XJj 



and general form of the wave function 

u 2 
u 3 

\ Ui ) 

the following system of equations will be obtained 

(E - m + jxB)ui - p z u 3 - (Il x - ill y )u 4 = 0, 

(E -m- fxB)u 2 + PzUa - (U x + iU y )u 3 = 0, 

(E + m - flB)u 3 - p z ui - (U x - itt y )u 2 = 0, 

(E + m + jxB)u 4 + p z u 2 — (H x + iRy)u\ = 0. 

To solve equations (1501) one can assume 

ui = CiI n (g), u 2 = C 2 I n -i(£,), U3 = C 3 I n (£), u 4 = C 4 I n -i(0i 
where the function is defined as 

I n (£) = N n e-? /2 H n (0; n 

Px 



(49) 



(50) 



(51) 



0,1,2,... 



f = VeBy + 



N„ 



eB 
14 



VeB 



(52) 



The function H n (£) is the Hermit polynomial and in according to the definition = 
0. Using this choice for normalization constant N n , it can be showed that the functions 
I n (£) satisfy the completeness relation 

oo 

^i n (0m') = 5(y-y'), (53) 

n=0 

where £' is given by replacing y by y' in the definition of the £ (1521) . After substitut- 
ing equations (IBTj) in the system of equations ( )50|) . the system of linear equations for 
coefficients Cj will be obtained as follows 

(E =f m ± A-B)Ci,3 - PzC 3 ,i - V2neBC 4)2 =0 , 

{E^m^ jlB)C 2A + Pz C i>2 - V2neBC 3A =0 . (54) 

Equations ( 154"|) have a non trivial solutions for Cj if the Determinant of the coefficients 
be zero. Imposing this condition yields the spectrum of the energy 

E n , s = yjp% + (V^ + 2neB - sfiB) 2 ] s = ±1. (55) 

In the energy spectrum ( |55l) index of n is called Landau level. To solve equations ( 1541) 
it is better to introduce the fermion spin projection operator on the direction of the 
magnetic field as 

fi z = rnL z + P{t xllV (56) 

Since the operator fi z commutes with the Dirac Hamiltonian Hp, the wave functions 
are eigenstate of this operator 

fizi) = fJ-oi>, (57) 

where /io is the eigenvalue of the operator fi z . Then using the equations ( 1571) and ( I5T1) . 
one leads to the following system of liner equations for the coefficients Ci 

(m =f /i )Ci j4 ± \/2neBC A + = 0, (m ± fi )C 2 , 3 ± V2neBC 3>2 = 0. (58) 

Again by imposing the condition that, Determinant of the coefficients must be zero, fio 
can be determined as 

Ho = sVm 2 + 2neB; s = ±1. (59) 

Now by using the equations (154 p . (15 8 p and the following normalization condition, the 
coefficients Cj can be determined 

£|Ci| 2 = l. (60) 
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Therefore the wave function can be obtained as 



U s (y,n,p$) 



N s e^U s (y,n,p^ 



V ( 



2 

2 in 2 Fzsj ± n—1 



±1. 



(61) 



(0 / 



where the quantities p zs and 7„ are defined as 

Pz 



Pz 



In 



V2neB 



+ Vm 2 + 2neB' 



m 



and the normalization constant N s of the wave function is 

1 



(62) 



(63) 



One can easily find the wave functions for the negative energy spinors by doing the 
similar procedure. In this case it is easier to determine the upper component with 
respect to the lower one. After some manipulations one has 

tfj s (x) = N s e ip - x *V s (y,n,w); 

( C-^Pzs - ^7n)4(£) \ 



(-^7n-^ S )4-l(0 
V C-^Pzsln + ^In-liO J 



£ = VeBy 



Pa 



r eB' 



(64) 



2 fzs in i 2 

where quantities p zs and 7 n and normalization constant N s have the same definition as 
the positive energy solutions and can be determined, respectively, from equations (I5"2l 
and fl63|) . Neutron has not charge but due to its anomalous magnetic moment its wave 
function and energy spectrum can be modified by the magnetic field. Solving Dirac 
equation for this case is simpler than proton case. By doing the same procedure as the 
proton case, the spectrum can be obtained as 



E s = \IP 2 Z + {ym 2 + p 2 ± - sfiB) 2 ; 
and the wave functions will be given by 

1),(x) = N s e- ip - x U s (p); 

I C-¥ + ^PzsPs) \ 



±1. 



(65) 



U s {p) 



(-^PzsPs + if 2 ) 
C-^Pzs + ^Ps) 
V C-PPs - ^Pzs) J 



±1, 



(66) 
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where quantities p zs , p s and p± are defined by 

Pz 



Pz 



Ps 



P X + SiPy, 



Ps 



Ps 



+ y?m 2 + p\ ' 



rn 



P± = VP+P- 



and the normalization constant given by 

1 



Pi 



'P+P- 



^{TTMsW+K) 

Finally the solutions for the negative energy case can be obtained as 

^ s {x) = N s e ip - x V s (p); 

( C-^Pzs + ^Ps) \ 

(I 

V S (P) = 



C-^Ps - ^Pzs) 
+ ^PzsPs) 

V (—■^rPzsPs + / 



(67) 



(6£ 



(69) 



2 fzsfs i 2 

Definition of the parameters in the wave functions (1691) are the same as the equations 
(EH and (EF 



B Average of the squared scattering matrix ele- 
ments 

The matrix elements that are sum of the neutral (1201) and the charge current amplitude 
(12TT) . contain the following integrals 

dye-Wv-^UOlmiU (70) 

where is given by replacing p x by p' x in the definition of the £ ( lo"2"|) . Details of this 
integration can be found in |l9j, however the final answer is as follows 



ZL e - J 7V2 e i^o ( + i% ) m - n L™- n ( V 2 ); m > n > 0, 
ml 



L n,m 




n > m > 0, 



(71) 



in which is the generalized Laguerre polynomial and variables rj and 0o ar e defined 

as 



Px-Px 

V2^B ' 



2 i 2 
Vx + Vy, 



k y ky 



(p'.+P^ik'y-ky) 

2eB 



(72) 
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By using the explicit form of the wave functions and the definition flTOj) . one can find 
the nonrelativistic leading order term in the amplitude as follows 



M = M nc + M cc , 
such that the neutral current amplitude matrix element is given by 

Mnc 



(73) 



n 2 



mt 



.1 + r. .1 + r\ ~, ~ .1 + r\ .1 — r\ ,~. ~ , 

c—)(-2-)(^ - kz) + (—)(-2-)(fc; - *;+) 



1 — r 1 -+- r ~ 1 — r 1 — r ~ 



x 



, .1 + s. .1 + s\ . . . .l + s.,1 — s' .1 — s w l + s' lT 

)(— 2~~ )(Pz -P*) + (— g— )(— )7n' - (— g— )(— )TnKn,n' 
{(— 2~ )(— J" )(& ~ ^~Y~^~T~^ n ' + 2~)7nKn-l,n'-l 



(74) 



and the charged current amplitude is given by 



2ZV 



, 1 + r 1 + s' r 1 + r . . 1 — s' . _ . 1 - r 1 + s' r 

„. — " — )W - *■) + ( — )( — )1W " ( — )( — )*-} 

7T 

r/ l + s w l + r' r, ~ x /I + 5 W 1 — r' ~ .1 — s w l + r\ „ . T 
x {(—)(— - P.) + ~ (— )(— )7»K»,n' 

, , 1 + r 1 — s' , r , 1 — r w 1 + s' . _, . 1 — r w 1 — s' s . . T N , 
+ {-(—)(-2-)fc+ + (—)(-2-)7n' - (—)(-2-)(p* " fc*)} 

r/ 1 + s. .1 + r' ~. . .l + s w l — r'.~ .1 — s w l+r' . T 

X {(—X—j-M -Pz) + - (— )(— )7»Kn,n'-l 

, 1 + r 1 + s' . j- 1 + r . , 1 - s' . _ , 1 — r . . 1 + s' . 7 , 
+ {(— )(-2~ )(P 2 - + (— )(-2~ )7n' - (— )(~2~ 

. ,l + s.,l — r'. , 1 — s. . 1 + r' ,~ .1 — s w l — r\,~. 

X {-( — )(— )7n + (— )(— J* 7 - " ( — )( — )(*; - P*)K-l,n' 
+ {(^)(^)* + " ( 1 T I )( i ^)7n' + (^)(^)(P1 - k)} 



(75) 



Using the following nonrelativistic expansion of the energy for proton as 



-En 



and for the neutron as 



m N + — ^— H s/i pJ B, 

imiv ^v 



£ s = m N + - s/j, n B, 

2m N 



(76) 



(77) 



along with the equations f )62|) and f )67|) lead to 

p z ~ V2neB ~ k z 

Pz — 7; j In — j kz 



2mjv 2m]y 2rriN 

with the similar definitions for the primed quantities. 
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